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ABSTRACT 


This paper presents a performance evaluation of two solutions to reduce 
computational burden of the traditional Weighted Least Squares Algorithm 
for power system state estimation: Simplified methods SWLSI / SWLS2 
based on full constant matrices and Fast decoupled FDWLS based on 
decoupled constant matrices. First, the algorithms were tested on IEEE 14 
and 118 bus transmission systems. Second, the solutions were tested on a 
rural distribution feeder to evaluate the response of the algorithms to high 
R/X ratio. Results show that for transmission systems, FDWLS is the fastest 
method but more sensitive to erroneous measurements. Simplifications 
considered in FDWLS, are not valid in distribution systems with high R/X 


Power system state estimation ratio this results in slowing down the algorithm convergence speed 


Weighted least squares considerably compared to SWLS2 which performs well. SWLS2 algorithm 
algorithm presents a promising solution to reduce computation time for application in 
future smart grid. 
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1. INTRODUCTION 

Recently, the power system is undergoing big changes due to the massive integration of renewable 
energy resources to transmission and distribution systems which leads to a complex bidirectional power flow. 
This conjucture incites researchers to improve the speed and reliability of state estimation algorithms to 
ensure an efficient real time monitoring of futur Smat Grid. 

The weighted least-squares (WLS) method is the most State Estimation Algorithm used in control 
centers all over the world. In fact, WLS state estimation algorithm provides the best estimation quality and 
good convergence rate. However, the gain and Jacobian matrices need to be recalculated each iteration which 
needs a large amount of calculation, a big memory requirement and long computing time [1]. One of 
the major solutions proposed in literature to circumvent the computational burden is the Fast-decoupled WLS 
(FDWLS) technique based on decoupled constant matrices [1-5] as used to speed up load flow 
calculation [6-8]. The fast-decoupled formulation has proven its efficiency to reduce computation time and 
data storage and has found wide acceptance in the industry, various versions have been implemented in 
control centers all over the world [4]. However, the decoupled method may fail to provide a solution on ill- 
conditioned systems as high R/X ratio of distribution branches or in presence of erroneous measurements. 
Therefore, it will be interesting to develop techniques which ensure a compromise between the reliability of 
full Weighted Least Squares Algorithm and Convergence speed of fast decoupled WLS and able to overcome 
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limitations discussed above. Authors of [5], proposed two variants of WLS method with full constant 
matrices (SWLS1 and SWLS2) based on Dishonest Gauss Newton Method used for Load Flow [2, 9-11] to 
reduce calculation time without decoupled simplifications. Though the principle of these methods is very 
popular, few studies analysed its convergence performances [5, 12-15] on Transmission systems and no study 
apply it before to a distribution system. This article is the first one exploring the performance of full versus 
decoupled constant matrices on distribution system state estimation with WLS algorithm. The rest of 
the paper is organized as follow: Section 2 presents a description of the traditional weighted least squares 
algorithm and the alternative solutions to reduce computational burden. The algorithms studied are tested 
both on transmission and distribution systems, employing accurate and erroneous measurements. Simulation 
results are presented and discussed in Section3. The main findings of the paper are as follows: 

- For Transmission Systems, FOWLS based on decoupled constant matrices is the best method, it is 2 to 
5 times faster than basic WLS and requires half storage capacity. 

- SWLS1I method considering constant Jacobian matrix only is not reliable and does not represents any 
advantage on reducing computation time. 

- For distribution systems, FDWLS based on decoupled constant matrices converges slowly and requires 
high iterations number while the SWLS2 algorithm evaluating the Jacobian and gain matrices once at 
the flat start is 2 times faster than basic WLS with the same reliability which make it less sensitive to 
erroneous measurements and high R/X ratio compared to FOWLS. Therefore, SWLS2 presents a good 
alternative to reduce computation time in future power systems. 


2. WEIGHTED LEAST SQUARES ALGORITHMS 
2.1. Basic algorithm 
The Network model employed is the single-phase model with N buses and m measurements 
gathered from remote meters. Most commonly used measurements are the line power flows, bus power 
injections and bus voltage magnitudes. 
The aim of state estimator is to provide the best possible values of the bus voltage magnitudes and 
angles by processing the available network data recognizing that there are errors in the measured quantities. 
The starting equation for the WLS state estimation algorithm is [4, 16-18]: 


z= h(x) +e (1) 


Where: z is the (mx1) measurement vector; 

x is an (nxl) state vector to be estimated: The number of estimated states is n=2*N-1, since 
the balance phase’s is already known O=0. 

e is an (mx1) measurement error vector. 

h is the vector of nonlinear functions that relate the states to the measurements defined: 

Real and Reactive power injection at bus 1: 


Pi = Vi die; Vj (GijcosOij + BijsinOij) (2) 
Qi = Vi diz; Vi(GijsinOij + Bijcos@ij) (3) 


Real and reactive power flow from bus 1 to bus j: 


Pij = Vi? (gsi + gij) — ViVj(gijcos6ij + bijsinðij) (4) 
Qij = —Vi? (bsi + bij) — ViVj(gijsinOij — bijcosdij) (5) 
Where: 


Vi is the voltage magnitude at bus 1 

ĝi is the phase angle at bus i 

dij = 0i — 0j 

Gij + jBij is the ijth element of the Y-bus matrix 

gij + jbij is the admittance of the series branch between bus i and bus J 

gsi + jbsj is the admittance of the shunt branch at bus 1. 

In practice, it is required to have the number of measurements larger than number of states, this is 
called redundancy [19]. So, state estimator can consider the various operation layouts used and to cover for 
the unavailability of transmission and telemetering equipment failures [20]. A measure of the redundancy 
may be denoted by the redundancy factor n, which 1s defined as [21]: 


Int J Pow Elec & Dri Syst, Vol. 11, No. 3, September 2020 : 1287 — 1297 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1289 


Dimension of Z m m 


Dimension of x n 2N-1 





The measurement errors ei are assumed to satisfy the following statistical properties: 
First, the errors have zero mean: E(ei) = 0, 1 = 1, ..., m. Second, the errors are assumed to be 
independent, such that the covariance matrix is diagonal: 


Cov(e) = E (e, e!) = R = diag {617, 627, ...,6m7} (7) 


The solution to the state estimation problem can be formulated as a minimization of the following 
objective function: 


JO = Er, RON fz- ney) Re — AG) 
(8) 


To find the minimization of this objective function the derivative should be set to zero. 
The derivative of the objective function is denoted by g(x): 


g(x) = Z2 = -H1 REZ — h(x)] = 0 (9) 


where: H(x)= Ch(x)/Ox, called the measurement Jacobian matrix. 
Ignoring the higher order terms of the Taylor series expansion of the derivative of the objective 
functions yields an iterative solution as shown below: 


x2 = xk + [GX] HH TRIT- h] (10) 


Where the gain matrix, G, is defined as: 
G(x") = “3 = H™R-1H (11) 


For the first iteration of the optimization problem, an initial guess has to be made for the state vector 
xO corresponding to the flat voltage profile, or flat start. A flat start refers to a state vector where all 
the voltage magnitudes are 1.0 per unit and all the voltage angles are O degrees. With respect to the state 
vector: the measurement functions, Jacobian and gain matrices are calculated each iteration until the absolute 
difference between two successive values of x is less than a chosen tolerance ¢. WLS flowchart presented in 
Figure 1, resumes the principal steps of the algorithm. 
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Figure 1. Basic WLS flowchart 
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2.2. Simplified algorithms with full constant matrices 

Mostly, it is observed that the elements of Jacobian and gain matrices do not significantly change 
between flat start initialization and the converged solution [5]. So, their evaluation could be restrained to 
some first iterations without altering the estimation quality. Based on this principle two methods have 
emerged: SWLS1 and SWLS2. 


2.2.1. The first simplified method SWLS1 

The first simplified method (SWLS1) calculates the Jacobian matrix at every iteration but preserves 
the gain matrix constant after an iteration k chosen as presented in Figure 2. The obvious advantage of this 
method is to reduce the number of calculations of the gain matrix. 


2.2.2. The second simplified method SWLS2 


The second method SWLS2 admits that the gain and the Jacobian matrices remain constant after an 
iteration k chosen as presented in Figure 3 [5, 13-15]. 
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Figure 2. SWLS1 flowchart Figure 3. SWLS2 flowchart 


2.3. Fast decoupled weighted least squares algorithm (FDWLS) 

The principle of Fast Decoupled State Estimation Techniques consists on exploiting 
the active/reactive decoupling property based on fast decoupled load flow methods [2, 9-11]. Indeed, for 
large scale power systems the transmission lines have a very high X/R ratio. In such a case, the real power 
changes are less sensitive to changes in voltage magnitude and can be ignored. Similarly, the reactive power 
change is less sensitive to changes in angles. Making these simplifications, the gain matrix in WLS state 
estimation algorithm can be simplified [10-11]. 

In the decoupled formulation, the measurement vector Z is portioned into two parts: active ZA and 
reactive ZR components [2-4]. 


o. Qinj 
Z=) Za h) Za (0r) a2 
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Where: 
Pinj, Qinj: are respectively real and reactive power injection measurements. 
Pflow, Qflow: are respectively real and reactive power flow measurements. 
V: Voltage measurements. 
ee the definition ae the : acobian matrix H, and covariance matrix R can be written as: 
H H 
[Hra Herrl’ Ra 
Based on P-O and Q-V decoupling, the off-diagonal blocks H4ę and Hp, in the measurement 
jacobian H are ignored: 


[Haa 0 
HO nl e 


Then the gain matrix is expressed as follow: 


G44 0 
“HD P (15) 
Where: 
Gaa = HaaRa Hag (16) 
Grr = HarRe Ure (17) 


Several versions of the fast-decoupled state estimator have been proposed in literature depending 
upon the assumptions and approximations adopted. The version recognized as presenting best performance 
has the following features [2]: 
The Jacobian and gain matrices are evaluated once, at the flat voltage start: all V=1 p.u and O = 0 
degrees. 

- The branch series resistances are ignored in forming the elements of the Jacobian H, which is equivalent 
to replacing branch susceptances bij by - 1 / xij. Where xij is the reactance of the series branch between 
bus i and bus J. 

- A transformed measurement vectors ZA’ and ZR’ are used by dividing the flow and injection 
measurements by the corresponding calculated voltage magnitude. 

- For better convergence characteristics, a block sequential solution scheme is used wherein: 


G,,A0 =TA (18) 
TA=HAAT*RA!*AZA ’ (19) 
AZA’=(ZA-hAV/V (20) 


hA is the vector of active measurements functions. 
A@ is solved and the updated 8's are used in the RHS of: 


GprpAV = TR (21) 
Which is then solved for AV. 

TR=HRR!*RR!*AZR’ (22) 
AZR’=(ZR-hR)/V (23) 


hR is the vector of reactive measurements functions. 
For more details, about how to calculate the elements of the fast-decoupled Jacobian, reader can 
refer to [2-4]. Algorithm steps are depicted in Figure 4: 
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Figure 4. FDWLS flowchart 


3. SIMULATION RESULTS 

The objective of this section is to compare the performance of the Weighted Least Squares 
Algorithms described in the previous Section: Basic WLS, SWLSI, SWLS2 and FDWLS in terms of 
computation time, convergence rate, estimation quality and number of iterations. 


3.1. Description of simulations 

The four algorithms have been tested on three study cases IEEE 14, 118 bus systems and a 12-bus 
rural distribution system under two scenarios: In scenario 1, accurate measurements equal to load flow results 
have been employed. In scenario 2, measurements were altered by a noise of 10% to evaluate the response of 
the studied algorithms to erroneous measurements. 

For IEEE bus systems, the network data files can be downloaded from Power Systems Test Case 
Archive [22]. 

To compare the state estimate accuracy of the following simulations, mean absolute percentage error 
(MAPE) is introduced as follow [23]: 


Xt—Xe 
Xt 


MAPE = = n x 100% (24) 








Where, Xt is the true value of system state obtained from load flow results and Xe is the estimated 
state. A smaller value of MAPE indicates a more accurate state estimation result. 

MAPEV: mean absolute percentage estimation error of voltage magnitude. 

MAPEO: mean absolute percentage estimation error of voltage angle. 

The simplified algorithms have been tested for different value of the iteration k from which just 
the gain matrix is considered constant SWLS1 or the gain and Jacobian matrices are constant (SWLS2). 
Convergence tolerance used in these tests is: e=10-4 for both voltage magnitude and voltage angle states. 


3.2. Test systems 

For all test cases, measurements were set to ensure a redundancy factor >1. They were chosen of 
different types and uniformly distributed through the network to ensure observability [16, 17]. Weight of all 
measurements is assumed 1. 


a. [EEE 14 bus system: 
For IEEE 14 bus system test case, a set of 41 measurements (n =1,5) is chosen as: 
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- 1 voltage magnitude at bus 1. 

- 8 real power injections and 8 reactive power injections at buses: 2, 3, 7, 8, 10, 11, 12, 14. 

- 12 real power flow and 12 reactive power flow on branches: 1-2, 2-3, 4-2, 4-7, 4-9, 5-2, 5-4, 5-6, 
6-13, 7-9, 11-6, 12-13. 


b. IEEE 118 bus system: 
A set of 726 measurements (7 =3) is chosen a follow: 
- 118 voltage magnitudes for all buses. 
- 118 real power injections and 118 reactive power injections for all buses. 
- 186 real power flow and 186 reactive power flow for all network branches. 


c. 12 bus distribution system: 

A physically existing rural distribution Feeder [24] is considered. Figure 5 shows the single-line 
diagram of the system, where node | is the substation. The system data including line data and load data can 
be found in [25]. 


S2 S3 S4 S5 S6 S7 Sg So S19 S11 S12 


Figure 5. Single-line diagram of the 12-bus distribution system. 


A set of 25 measurements (n =1,1) is chosen as follow: 
- 1 voltage magnitude at bus 1. 
- 1 real power injection and 1 reactive power injection at bus 1. 
- 11 real power flow and 11 reactive power flow for all network branches. 


3.3. Results of scenario 1 (accurate measurements) 
In this section, WLS algorithms have been tested with accurate measurements equal to load 
flow results. 


3.3.1. Simulation results for IEEE 14 bus system 

As seen in Table 1: 

FDWLS is the fastest method even if number of iterations is increased compared to the basic WLS. 
FDWLS solution changes slightly but precision remains good. 

SWLS1I does not converge at flat start (k=1), for k=2 number of iterations and computation time is 
increased compared to the other methods. The solution is the same as Basic WLS. 

SWLS2 applied at flat start (k=1) converges in half time compared to the basic WLS with a good 
quality of estimation even if number of iterations required is important compared to the other methods. 


Table 1. Performance evaluation of WLS algorithms for IEEE 14 bus system under scenariol 


Algorithm Computation time (seconds) Iterations number MAPEV (%) MAPEO (%) 
Basic WLS 0.004541 4 1.28 2.41 
SWLS1 (k=1) Program does not converge 
SWLS1 (k=2) 0.005751 7 1.28 2.41 
SWLS1 (k=3) 0.004102 4 1.28 2.41 
SWLS2 (k=1) 0.002113 8 0.57 0.87 
SWLS2 (k=2) 0.002550 4 1.33 ZS 
SWLS2 (k=3) 0.003545 4 1.28 2.41 
FDWLS 0.001418 6 0.34 0.45 


Figure 6 shows that FOWLS and SWLS2 (k=1) are closer to the true value than basic WLS, 
the same figure appearance was obtained for voltage angle. 
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Figure 6. Estimated voltage magnitude (P.U) by different WLS algorithms for IEEE 14 bus 
system under scenario 1 


3.3.2. Simulation results for IEEE 118 bus system 

According to Table 2: here also, SWLS1 does not converge for k=1. It seems that constant gain 
matrix associated to variable Jacobian may lead to convergence problems. FOWLS computation time is the 
lowest. The all methods have provided a good estimation quality. 


Table 2. Performance evaluation of WLS algorithms for IEEE 118 bus system under scenariol 


Algorithm Computation time (seconds) Iterations number MAPEV (%) MAPEO(%) 

Basic WLS 0.182347 4 1.79 1.74 
SWLS1 (k=1) Program does not converge 
SWLS1 (k=2) 0.265629 7 1.79 1.74 
SWLS1 (k=3) 0.156127 4 1.79 1.74 
SWLS2 (k=1) 0.084726 10 1.81 1.13 
SWLS2 (k=2) 0.098518 5 1.77 1.72 
SWLS2 (k=3) 0.135495 4 1.79 1.74 

FDWLS 0.037558 5 2.07 1.1 


3.3.3. Simulation results for 12 bus distribution system 

As noticed in Table 3, SWLS1 has converged for k=1 because the final state is too close to 
the initial state (Basic WLS converges only in 3 iterations). 

FDWLS convergence speed slowed down considerably compared to the other algorithms. 
The solution of FOWLS is the less accurate, and the program required a high iteration number (60,5 against 3 
for the basic formulation). SWLS2 (k=1), converged in less time and its solution is accurate as the basic 
WLS ones. 


Table 3. Performance evaluation of WLS algorithms for 12 bus distribution system under scenario 1 


Algorithm Computation time (seconds) Iterations number MAPEV (%) MAPEO (%) 
Basic WLS 0.002201 3 0.01 1.02 
SWLS1 (k=1) 0.003601 5 0.01 1.02 
SWLS1 (k=2) 0.002021 3 0.01 1.02 
SWLS2 (k=1) 0.001545 5 0.01 1.02 
SWLS2 (k=2) 0.001906 3 0.01 1.02 
FDWLS 0.003879 60.5 0.76 19.87 


Results are depicted in Figure 7: the differences between estimated voltages by SWLS2(k=1), Basic 
WLS and the true values are so small that they cannot be distinguished. 
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Figure 7. Estimated voltage magnitude (P.U) by different WLS algorithms for 12 bus distribution system 
under scenario 1 


3.4. Results of scenario 2 (erroneous measurements) 
In this section, measurements were altered by a noise of 10% to evaluate the response of the studied 
algorithms to erroneous measurements. 


3.4.1. Simulation results for IEEE 14 bus system 

Table 4 shows that bad measurements alter the estimation quality for all WLS algorithms. Specially, 
FDWLS which provided the less accurate solution. However, FOWLS method still constantly the fastest one. 
For this case, SWLS1 have failed not only at flat start but also for the second iteration (k=2). 

According to Figure 8, All WLS algorithms solutions are far from the true values. The FDWLS 
curve is the furthest while Basic WLS curve is the nearest one to the true values curve. Those results are 
the inverse of those obtained in scenario 1. 


Table 4. Performance evaluation of WLS Algorithms for IEEE 14 bus system under scenario2 


Algorithm Computation time (seconds) Iterations number MAPEV (%) MAPEO (%) 
Basic WLS 0.005858 5 7.82 5.07 
SWLS1 (k=1) Program does not converge 
SWLS1 (k=2) Program does not converge 
SWLS1 (k=3) 0.005266 5 7.82 5.07 
SWLS2 (k=1) 0.002391 8 9.08 7.32 
SWLS2 (k=2) 0.003127 5 7.3 4.17 
SWLS2 (k=3) 0.003545 4 7.8 5.04 
FDWLS 0.001459 6.5 9.6 8.15 





Figure 8. Estimated voltage magnitude (p.u) by different WLS algorithms for IEEE 14 bus system under 
scenario2 


3.4.2. Simulation results for IEEE 118 bus system 

Results presented in Table 5, are less affected by noise compared to the previous case IEEE 14 bus 
system, especially for voltage angle state. In fact, the set of measurements chosen for this studied model, 
presents a high redundancy level (n =3) which allowed to the WLS Algorithms to provide good estimation. 
SWLS1 has converged for k=2 in more iterations and computing time compared to the other algorithms. 
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Table 5. Performance Evaluation of WLS Algorithms for IEEE 118 Bus system under Scenario2 


Algorithm Computation time (seconds) Iterations number MAPEV (%) MAPEO (%) 
Basic WLS 0.184747 4 6.65 0.99 
SWLSI1 (k=1) Program does not converge 
SWLS1 (k=2) 0.456256 14 6.66 0.99 
SWLS1 (k=3) 0.170388 4 6.65 0.99 
SWLS2 (k=1) 0.084137 10 7.21 1.85 
SWLS2 (k=2) 0.095147 5 6.54 0.87 
SWLS2 (k=3) 0.130684 4 6.66 0.99 
FDWLS 0.041343 5.5 8.37 2.93 


3.4.3. Simulation results for 12 bus distribution system 
As noticed in Table 6, here also the FDWLS has converged slowly with high iteration numbers. 
While, SWLS2(k=1), converges in reduced time compared to the basic WLS, with the same reliability. 


Table 6. Performance evaluation of WLS algorithms for 12 bus distribution system under scenario2 


Algorithm Computation time (seconds) Iterations number MAPEV (%) MAPEO (%) 
Basic WLS 0.002934 4 10.35 9.07 
SWLS1 (k=1) 0.005536 9 10.35 9.06 
SWLS1 (k=2) 0.002801 4 10.35 9.07 
SWLS2 (k=1) 0.001602 6 10.34 9.08 
SWLS2 (k=2) 0.002208 4 10.35 9.07 
FDWLS 0.004408 68.5 9.6 9.11 


4. DISCUSSION 

Applied on transmission systems with low R/X ratio, FOWLS is the fastest WLS Algorithm. It is 2- 
5 times faster than basic WLS even if convergence requires half to two iterations more than basic WLS 
(simplifications assumed diminish the true quadratic convergence properties of basic WLS). FDWLS 
requires half storage capacity since off diagonal elements of H and G matrices are neglected. Diagonal 
Jacobian elements are simples and contain only reactance branches. However, applied on distribution system 
with high R/X ratio, FOWLS was the slowest algorithm and required a huge iterations number because 
assumptions considered are based on transmission system features (active/reactive decoupling). 

SWLS1 algorithm does not converge for first iterations. It seems that constant gain matrix 
associated to variable Jacobian may lead to convergence problems. SWLS1 is not reliable and does not 
present any advantages on reducing computation time. 

SWLS2 algorithm applied at first iteration (k=1), requires the evaluation of gain and Jacobian 
matrices just once which reduces computation time. Indeed, SWLS2 (k=1) is 2 times faster than basic WLS 
even it requires a higher iterations number. SWLS2(k=1) has the same characteristics of basic WLS which 
make it less sensitive to erroneous measurements and high R/X ratio compared to FOWLS and it could be 
applied in distribution network state estimation. 


5. CONCLUSION 

This paper has presented a performance comparison between two solutions to reduce computational 
burden of traditional Weighted Least Squares (WLS) Algorithm: Solution 1 based on full constant matrices 
(simplified methods SWLSI / SWLS2) and Solution 2 based on decoupled constant matrices (Fast 
Decoupled WLS FDWLS). Algorithms have been tested on three study cases IEEE 14, 118 bus systems and 
a 12-bus rural distribution system. For each case, simulations were performed on accurate and erroneous 
measurements. 

Simulation results have shown that the performance of FDWLS algorithm depends on the 
application environment: applied on transmission systems, FDWLS is the fastest WLS Algorithm. However, 
on distribution system with high R/X ratio, FOWLS was the slowest algorithm and required a huge iterations 
number. Therfore, FDWLS is not suitable for distribution system state estimation. 

SWLS2 (k=1) is 2 times faster than basic WLS with the same reliability, which make it less 
sensitive to erroneous measurements and high R/X ratio compared to FOWLS. Therefore, SWLS2 applied at 
the flat start presents a good alternative to reduce computation time in future distribution system 
state estimation. 
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